Draft version September 15, 2011 

Preprint typeset using I^'T^]X style emulatcapj v. 11/10/09 



MODELING THE EFFECTS OF STAR FORMATION HISTORIES ON Ha AND ULTRA-VIOLET FLUXES IN 

NEARBY DWARF GALAXIES 

Daniel R. Weisz^'^, Benjamin D. Johnson^, L. Clifton Johnson^, Evan D. Skillman^, Janice C. Lee* ''', Robert C. 
Kennicutt, Jr.^, Daniela Calzetti®, Liese van Zee^, Matthew S. Bothwell"^, Julianne J. Dalcanton\ Daniel A. 

Dale*^, Benjamin F. Williams"' 

Draft version September 15, 2011 



o 

(N 

CD 

m 



O 

U 

6 



> 

in 
o 

a^ 

(N 

0^ 
O 



ABSTRACT 

We consider the effects of non-constant star formation histories (SFHs) on Ha and GALEX far 
ultra-violet (FUV) star formation rate (SFR) indicators. Under the assumption of a fully populated 
Chabrier IMF, we compare the distribution of Ha-to-FUV flux ratios from ~ 1500 simple, periodic 
model SFHs with observations of 185 galaxies from the Spitzer Local Volume Legacy survey. We find 
a set of SFH models that are well matched to the data, such that more massive galaxies are best 
characterized by nearly constant SFHs, while low mass systems experience bursts amplitudes of 
30 (i.e., an increase in the SFR by a factor of 30 over the SFR during the inter-burst period), burst 
durations of tens of Myr, and periods of ~ 250 Myr; these SFHs are broadly consistent with the 
increased stochastic star formation expected in systems with lower SFRs. We analyze the predicted 
temporal evolution of galaxy stellar mass, i?-band surface brightness, Ha-derived SFR, and blue 
luminosity, and find that they provide a reasonable match to observed flux distributions. We find that 
our model SFHs are generally able to reproduce both the observed systematic decline and increased 
scatter in Ha-to-FUV ratios toward low mass systems, without invoking other physical mechanisms. 
We also compare our predictions with those from the Integrated Galactic IMF theory with a constant 
SFR. We find that while both predict a systematic decline in the observed ratios, only the time variable 
SFH models are capable of producing the observed population of low mass galaxies (Af* < 10^ M©) 
with normal Ha-to-FUV ratios. These results demonstrate that a variable IMF alone has difficulty 
explaining the observed scatter in the Ha-to-FUV ratios. We conclude by considering the limitations 
of the model SFHs, and discuss the use of additional empirical constraints to improve future SFH 
modeling efforts. 

Subject headings: galaxies: dwarf 
formation 



galaxies: evolution — galaxies: formation — galaxies: star 



1. INTRODUCTION 



Two of the most widely used tracers of recent star for- 
mation arc the nebular Ha recombination line and the 
ultra-violet (UV) continuum. The Ha arises from recom- 
bination of gas ionized by photons from massive stars (> 
15 Mq) and is expected to be observed over the typical 
lifetimes of extremely massive stars (< 5 Myr). The UV 
continuum is due to non-ionizing photospheric emission 
from stars with M > 3 Mq, which have lifetimes < 300 
Myr. In tandem, these two integrated tracers provide 
leverage on current and recent star formation in both 
nearby and distant galaxies (see lKennicutti[l998l and ref- 
erences therein). 

In principle, observed Ha and UV luminosities should 
yield consistent measures of star formation. Under the 
assumption of a constant SFH over a sufficiently long 
timeline (e.g., ^ 1 Gyr), the expected ratio of Ha and 
UV star formation rates (SFRs) should be co nstant with 
respect to both time and environment (see iKennicufrB 
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119981 and references therein) , assuming a fully populated 
and universal IMF, solar metallicity, all Lyman contin- 
uum photons ionize hydrogen, and with no attenuation 
due to dust. Deviations from this fiducial ratio would 
suggest that one or more of the underlying assumptions 
are not correct. 

A number of studies have demonstrated discrepancies 
in SFRs as measured by Ha and UV luminosities ( e.g. , 
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Iglesias-Paramo et al.l 120041 ). Recently, several studies 
have unambiguously shown that the observed discrep- 
ancy in Ha and UV SFR indicators is systematic, such 
that the observed ratio of Ha-to-U V flux declines with 
decreasing galaxy luminosity (e.g., iMeurer et al.l 120091 : 
iLee et al.ll2009bl : 'lBoseni et al.|[2009[ ) 



Despite extensive research, there remains no consen- 
sus for the cause of this trend. Independent studies 
have verified that factors such as assumed metallicity and 
choice of stellar m odels cannot be responsible for the ob- 
served trend (e^ lMeurer et all 120091 : ILee et al.|[2009bl : 
iBoselli et al.l |2009[ ). while investigations of other con- 
tributing factors, including non-constant SFHs, a vari- 
able IMF, and ionizing photon leakage have yet to yield 
conclusive results. 

In this paper, we undertake a focused investigation on 
the impact of non-constant SFHs on observed Ha-to-UV 
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Table 1 

Basic Properties of the LVL Galaxies 



Galaxy 


log{A/*) 


log(Ejj) 


Mb 


log[SFR(Ha)\ 


, F{Ha) 

iOg F{FUV) 


- K Source 


Name 


log{M0) 


log{L0 fcpc-2) 




(Mq yr-i) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


UGCA292 


5.65 


6.91 


-11.42 


-2.755 


0.031 


1 


UGC5364 


5.82 


7.24 


-11.40 


-4.135 


-0.739 


1 


UGC8091 


5.84 


7.54 


-12.07 


-2.700 


-0.007 


1 


UGCA438 


6.28 


7.65 


-12.35 


-4.316 


-1.668 


2 


UGC8833 


6.44 


7.60 


-12.42 


-3.331 


-0.468 


1 


UGC4483 


6.50 


7.72 


-12.86 


-2.446 


0.016 


2 


UGC9128 


6.51 


7.56 


-12.39 


-3.970 


-0.948 


1 


CGCG269-049 


6.52 


7.88 


-12.35 


-3.120 


-0.290 


1 


KKH37 


6.53 


7.56 


-11.67 


-3.818 


-0.319 


2 


UGCA281 


6.54 


8.32 


-13.50 


-1.391 


0.369 


1 



Note. — The observed and derived properties of the 185 galaxies considered in this study. For 
brevity, we list 10 here and make the rest available in machine readable format. The measurement 
of galaxy stellar masses and i?-band surface brightnesses are detailed in [[2] The H« SFR (column 
(4) have not been corrected for att enuation, while the Ha-to-FUV ratios in column (6) have been 
dust corrected as detailed in i|2.1l The values in this column are listed relative to the value of 
fiducial, whose value is reflected by the constant k = —13.17. Column (7) indicate s the source of 
the optic al observations, either from SDSS DR7 l lAbazaiian et ani2009l ) (1) or LVL llvan Zee et al.l 
lin prep.1) (2). The full data table will be available with the journal version of the paper or upon 
request. 

ratios. Our main objective is to find a set of simple SFH 
models (i.e., periodic bursts of star formation superim- 
posed on a baseline constant SFR) that are well matched 
to the observed distribution of Ha-to-UV SFR ratios, 
while being consistent with other available data. To do 
this, we compare the observed distribution of Ha-to-FUV 
ratios to the predicted distributions from a set of simple 
SFH models, in bins of galaxy stellar mass. Using the 
results of this comparison, we examine the ability of the 
highest probability model SFHs to explain the observed 
trend in Ha-to-FUV flux ratios versus stellar mass, R- 
band surface brightness, Ha luminosity, and Mb- In 
addition, we compare SFH model parameters and pre- 
dictions with previous studies that have analyzed trends 
in Ha-to-UV ratios using different SFHs or IMF assump- 
tions. We conclude the paper by discussing the strengths 
and limitations of simple model SFHs and suggest spe- 
cific ways to improve future SFH modeling efforts. 

2. THE OBSERVATIONAL DATA 

We test our SFH models on a sample of 185 nearby 
star-forming galaxies, whose observable properties are 
listed in Table [T] To arrive at this particular sam- 
ple, we began by considering the 390 galaxies with Ha 
and GALEX FTJV observations from the 11 Mpc Ha 
and UV Galaxy Survey CllHUGS: lKennicutt et alMOO^ 
iLee et ahl 1201 1[ ) . We further restricted the sample to 
those galaxies that also had high quality ancillary data 
to To allow us correct for the effects of dust and study 
the Ha-to-FUV ratio trends versus galaxy stellar mass 
and surface brightness. Specifically, we selected IIHUGS 
galaxies that are also members of the s maller Spitzer Lo - 
cal Volume Legacy survey (LVL; e.g.. lDale et all 120091) . 
The LVL sample contains 258 galaxies with comprehen- 
sive Spitzer IR imaging. However, only 59 galaxies in the 
sam ple also have homogene ous optical broadband imag- 
ing ( van Zee et aLllTn prep.t ). To increase the sample size, 
we also considered L VL galaxies that fall wit hin the foot- 
print of SDSS DR7 HAbazaiian et al.ll2009D . There are 
126 additional galaxies that met this requirement. Thus, 



we have a final sample of 185 galaxies with observations 
of Ha and GALEX FUV from IIHUGS, Spitzer IR from 
LVL, and optical broadband imaging from either SDSS 
or LVL. For consistency with the LVL optical imaging. 
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Fi gure 1. Ha-to-FU V ratios plotted versus the Ha SFR and Mb for 
the I Lee et al.| l |2009b|) sample (grey and navy points) and the sample 
considered in this study (navy points only). All ratios have been cor- 
rected for both foreground and internal extinction using the method of 
ILee et ah (2009b). To help illustrate the differences in the flux ratios, 
a constant k has been added. The value of k is —13.17, the negative 
of the expected flux ratios from our models for a constant SFH. The 
black dot-dashed line indicates the expected fiducial value for a con- 
stant SFH, solar metallicity, and fully populated Chabrier IMF. The 
two samples are generally consistent, although there is a noticeable 
deficiency in luminous galaxies {Mb ^ —16) with Ha-to-FUV ratios 
above the fiducial in the present study, relative to the larger sample of 
ILee et al.| | |2009bD . 
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Table 2 

Ha-FUV Flux Ratio Distribution Statistics for Stellar Mass Bins 



log SFRiHa) 
(1) 



SFR(Ha) 
SFR{FUV) 
(2) 



l(T-scatter 
(3) 



(4) 



log 



SFR(Ha) 
SFR(FUV) 
(5) 



l(T-scatter 
(6) 



loe 



SFR(Ha) 
SFR(FUV) 
(7) 



l(T-scatter 
(8) 



(9) 



0.50 


-0.14 


0.29 


20 


-0.32 


0.16 


-0.06 


0.05 


7 


-0.25 


-0.09 


0.20 


26 


-0.28 


0.13 


-0.03 


0.08 


22 


-0.75 


-0.13 


0.22 


44 


-0.22 


0.15 


-0.06 


0.10 


21 


-1.25 


-0.13 


0.22 


51 


-0.19 


0.15 


-0.06 


0.14 


36 


-1.75 


-0.18 


0.17 


58 


-0.24 


0.26 


-0.14 


0.24 


38 


-2.25 


-0.27 


0.22 


56 


-0.38 


0.45 


-0.26 


0.44 


40 


-2.75 


-0.51 


0.25 


31 


-0.38 


0.45 


-0.26 


0.44 


15 


-3.50 


-0.55 


0.57 


21 


-0.69 


0.31 


-0.53 


0.26 


6 



Note. — A c omparison between the Ha-to-FUV ratios as a function of Ha SFR between the current sample and that of 
I Lee et al.| l |2009b|) . Columns (l)-(4) are the du s t corre cted Hct SFRs, Ha-to-FUV ratios, Icr-scatter in the ratio and number of 
galaxies per bin taken from Table 2 in lLee et al.| | |2009bD . Columns (5) and (6) are the Ha-to-FUV ratio and Icr scatter for galaxies 
in the present study, with the extinction corrections of lLee et al.| l |2009b|) applied. Columns (7) and (8) are the Ha-to-FUV ratio 
and Icr scatter for galaxies in the present study, with the extinction corrections from fKennic utt et all 120091) and lHao etaT] 120111) 
applied. Column (9) is the number of galaxies per bin for the sample of 185 galaxies we consider in this study. 



fluxes measured from SDSS images were converted to the 
Johnson filter system using the transformations provided 
in Blanton ct al. ( 2005). 

In Figure [U we s how a com parison between the Ha- 
to-FUV ratios from iLee et all ((2009b' ) (navy and grey 
points) with those of the 185 galaxies considered in this 
study (navy points only). The Ha, FUV, and B-band 
fluxes, as well as the adopted dust corr e ctions in this 
figure have all been taken from ILee et alj (|2009bf ). 

Although similar in most aspects, there are two no- 
table difference between the samples. First, our subset 
of 185 galaxies are noticeably deficient in luminous sys- 
tems {Mb ^ —16) with Ha-to-FUV ratios significantly 
greater than the fiducial. Second, our sample contains 
fewer low luminosity systems than the larger sample of 
ILee et all (|2009b[ ). In both cases, the limiting factor is 
the lack of sufficient ancillary data. In the case of the low 
luminosity galaxies, nearly all that have been excluded 
from the present study lack sufficiently homogenous op- 
tical broadband coverage. For more luminous galaxies, 
some lack either ancillary optical imaging or fall out- 
side the boundaries of the LVL survey and therefore lack 
Spitzer imaging. 

We next examine the differences in Ha-to-FUV ra- 
tios between the two samples. In Table [21 we quantify 
the Ha-to-FUV ratios as a function of Ha SFR, after 
applying identical extinction corrections from ILee et ahl 
(j2009b() (see §??). As expected, the largest difference 
is seen at high Ha SFR, where our sample lacks galax- 
ies with high Ha-to-FUV ratios. For log[SFR{Ha)] > 
— 1.75, the mean Ha-to-FUV ratios in the present study 
are sy stematical l y ^0.2 dex lower than those in the sam- 
ple of ILee et al.l (|2009bO . Similarly, the la scatter over 
the same range is ~ 0.1 dex lower in the present sam- 
ple. However, for galaxies with \og[SFR{Ha)] < —1.75, 
the Ha-to-FUV ratios have comparable mean values and 
scatter, per bin of Ha SFR. Outside of the differences 
between the samples at high Ha SFRs, the two samples 
are in general agreement. This indicates that possible se- 
lection effects associated with the smaller sample in this 
study are minimal. 

2.1. Adopted Extinction Corrections 

For accurate comparison of Ha and FUV SFR indica- 
tors, corrections are required for extinction due to both 



Galactic foreground and to internal effects in each galaxy. 
We c orrect for foreg round extinction using the maps of 
Schle geTet alJ (119981). ass u ming the Milky Way extinc- 
tion law of ICardelli et al.l (11989). For internal attenua- 
tion, we use two types of corrections: the first based on 
direct IR measurement of reprocessed light and another 
usin g correlation- based prescriptions (e.g., see Appendix 
D in lLerov et alll200l . 

For ^ 70% of the sample, we have estimates of the 
total- IR flux based on Spitzer/MIPS far-IR (FIR) mea- 
surements, allowing us to use the 'ene rgy balance' ex- 
tinctio n correction m ethods derived by iKennicutt et al.] 
([200l for Ha and bv lHao et all (IMlh for FUV. These 
methods derive an intrinsic luminosity from the com- 
bination of unobscured (e.g.. Ha, FUV) and obscured 
(i.e., FIR continuum) signatures of star formation. From 
the intrinsic luminosity one can correct for the effects 
of internal dust attenuation on the observed Ha and 
FUV fluxes. For the - 30% of the sample where FIR 
measurements were not available, we used the same 
correlation-based attenuation corrections as iLee et aP 
(2()09b), namely a B-band luminosity prescription for 
Ana and a dust-law prescription for Apuv- While both 
methods are anchored to Balmer decrement extinction 
corrections, the self-consistency of the 'energy balance' 
method along with the availability of IR for the majority 
of our sample make it an appealing choice for this study. 

More explicitly, the 'energy balance' ext inction correc- 
tion for Ana from IKennicutt et al.l (|2009f) is: 

^ A 0.0019 L(TIR)\ 
Ag^ = 2.51og 1-f I ^ 1 

and from iHao et al.l (|2011l ). the correction for the FUV 
is given as: 



A 



FUV 



= 2.5 log 1 



0.37 L(TIR)\ 
L{FUV)obs ) 



(2) 



where the TIR luminosity is a MlPS-based measur ement 
made using the conversion in lDale fc Heloul (|2002[ ): 



L(TIR) = 1.559i//^(24/im) + 0.7686j//,.(70yLtTO) 

-|-1.347j//^(160/ito) 



(3) 



Overall, our internal extinction corrections yield simi- 
lar results to those used bv iLee et al.i (|2009bl . As shown 
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Figure 2. Same as F igure [T] except that we hav e applied the dust 
correction of Ke nni c utt et aLI II2009I ') for Ha and lHao et al.l ||2011I '| 
for the FUV. We adopt these dust corrections for all subsequent 
analysis in this paper. 



in Table [2l the main difference between the two sets of 
corrections is s een for the most himinous galaxies in the 
sample. Here, iLee et all ([2009b) find the mean Ha-to- 
FUV ratio to be ~ 0.2-0.3 dex below the expected fidu- 
cial, whereas the energy balance dust corrections place 
the mean ratios within ^ 0.05 dex of the fiducial. How- 
ever, for galaxies with low Ha SFRs, the adopted extinc- 
tion correction does not drastically influence the mean or 
scatter in the Ha-to-FUV ratios; these galaxies have low 
dust contents, such that extinction corrections have only 
modest impacts on their luminosities. We plot the flux 
ratios for our sample with the energy balance extinction 
corrections applied in Figure [2l 

Although differences in dust corrections are both im- 
portant and interesting in their own right, analyzing the 
impact of various extinction corrections on SFR indica- 
tors is beyond the scope of the present analys is . Instead , 
we refer the reader to IL. C. Johnson et all din prep.D . 
which investigates the reliability and consistency of vari- 
ous H a and FU V extinc ti on corrections, in c luding those 
from iLee et al.l (l2009bD. IKennicutt et al.l (|2009t) . and 
IHao et al.l (12'oTlh . 

2.2. Deriving Galaxy Masses 

The stellar mass for each galaxy in the sample was 
derived by spectral energy distributio n (SED) fitting as 
detailed in IB. Johnson et ahl fin prep.'), which closely fol- 
lows the methods detailed in Salini et al. (2007). Briefly, 
we fit the observed UV, optical data, and IR luminosities 
of each galaxy with a suite of Bruzual &: CharlotI (|2003D 
population synthesis models, whose parameters span a 
range in stellar metallicity, age, exponentially declining 
SFHs, and dust properties. For each model, we derive 
a stellar mass from the best fitting normalization of the 
model SED. We also construct the cumulative distribu- 
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Figure 3. Same as Figure [T] only plotted versus galaxy stel- 
lar mass and ij— band surface brightnesses (S^). In the lower 
panel, we have over -plotted the line fit line to the data studied by 
IMeurer et"an II2009I ') as the red dotted line. The Ha-to-FUV ratios 
considered in this study do not show as strong a correlation with 



tion function for the stellar mass from the likelihoods of 
the normalized model SEDs for each sample galaxy. The 
stellar masses reported here are the medians of the cumu- 
lative distribution function. In practice, the stellar mass 
is surprisingly robust, as the effects of dust and age on 
the optical mass-to-light ratio have corre spondingly sim- 
ilar e ffects on the observed optical color ([Bell &: de Jond 
120011 ). while the effects of stellar metallicity on the opti- 
cal mass-to-light ratios are minimal. The derived stellar 
masses for each system are presented in Table [T] 

The Ha-to-FUV ratio as a function of galaxy stellar 
mass (Figure [3]) shows subtle, but important differences 
compared to when the ratios are plotted versus Ha SFR 
or Mb, as in Figured For example, when considering the 
ratios versus Ha SFR, we see there are no galaxies with 
Ha-to-FUV ratios above the fiducial for \og[Ha{SFR)] 

< —2.5. Additionally, there are only a few galaxies with 
Ha-to-FUV ratios near or above the fiducial for systems 
with Mb ^ —13. From these two panels, one could con- 
clude that galaxies with low luminosities or low current 
SFRs must have systematically low Ha-to-FUV ratios. 

However, when plotting the same ratios versus stellar 
mass, one could reach a different conclusion. As shown in 
the top panel of Figure [31 a number galaxies with masses 

< 10^ Mq have Ha-to-FUV ratios near or above the fidu- 
cial. Recent bursts of SF have little impact on a galaxy's 
total mass, but can strongl y influence the obser ved Ha 
or optical luminosities (e.g.. iBell fc de Joii5l200lD . Thus, 
considering the data versus total stellar mass minimizes 
the covariances that affect the Ha-to-FUV ratios when 
plotted versus luminosity. 

2.3. Measuring R-band Surface Brightness 
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Table 3 

Ha-FUV Flux Ratio Distribution Statistics 
for Stellar Mass Bins 



Bin 
log(Af©) 
(1) 



, SFRjHa) 
°§ SFR(FUV) 
Median 
(2) 



(3) 



log(M^) >10.0 


-U.UO_Q Qg 


16 


10.0>Iog(M*)>9.0 




20 


9.0>Iog(M*)>8.0 




55 


8.0>Iog(M*)>7.0 


'^■^'-0.34 


71 


log(M*) <7.0 


r, 97+0-36 
-u.z/ _Q 32 


23 



Note. — The median Ha-to-FUV ratio per 
stellar mass bin. The errors bars correspond 
to the 16th and 84th percentile values of the 
distribution of observed ratios. The number of 
galaxies in each bin is listed in column (3). 



A similar study conducted by IMeurer et al.l (|2OO90 
found a strong trend between i?-band surface bright- 
ness (Sij) and Ha-to-FUV ratios. To facilitate a direct 
comparison, we utilize i?-band imaging from either SDSS 
or LVL, and follow the surface brigh t ness m easurement 
procedure described in IMeurer et al.l (j2006f) . Below, we 
briefly summarize the methodology. 

We first masked foreground stars and background 
galaxies identic ally to the method described in 
iDale et all ()2009[ ). With these objects removed, we com- 
puted the surface brightness profiles using elliptical aper- 
tures. The major and minor axes of a measurement 
aperture are scaled, preserving ellipticity, such that the 
ellipse encloses 50% of the total fiux. We then normal- 
ized the total fiux to the area of the aperture to de- 
rive the surface brightness. Particular care was taken 
to ensure that these scaled measurement apertures, con- 
structed using ellipticity and position data compiled from 
NED/SIMBAD, were properly s haped and centered. For 
consistency with IMeurer et al.l ()2009l ) , we note that the 
fiux percentage choice of 50% rather than 90% does not 
affect the overall trends in our data, beyond shifting to 
slightly higher surface brightnesses. 

We plot the Ha-to-FUV ratios versus i?-band surface 
brightness (Sfl) in the bottom panel of Figure [H Again, 
we see a decline in Ha-to-FUV ratios toward lower sur- 
face brightness systems. In this case, the systematic de- 
cline of Ha-to-FUV ratios is not as strong as when the 
sample is plotted versus Ha SFR, but is still more evi- 
dent than when the ratios are considered versus galaxy 
stellar mass. Like all optical luminosities, i?-band fiuxes 
are known to vary d ue to episodes of recent SF (e.g., 
iBell &: de JonEl[2001f) , which likely explains at least some 
of the correlation between Ha-to-FUV ratios and E^,. 
We note that the observed trend between Ha-to-FUV 
ratio and i?-band surface brightnes s in our sample is no t 
quite as strong as that presented in IMeurer et al.l (|2009l ). 
which we have over-plotted as the red dotted line in Fig- 
ure 131 We discuss some of the potential reasons for this 
disagreement in 



3. MODEL STAR FORMATION HISTORIES 

The model SFHs we consider in this analysis are sim- 
ple toy models that are primarily composed of periodic 
bursts superimposed onto a constant SFR. Given that 



modest variations in the period, duration, and ampli- 
tude of bursts can result in significantly different Ha, 
UV, and broadband fluxes, we consider a diverse and ex- 
tensive set of ~ 1500 model SFHs, whose properties and 
construction we describe below. 

3.1. Construction of the Model Star Formation 
Histories and Associated Fluxes 

We model SFHs that focus on two separate regimes: 
the ancient epoch (> 1 Gyr ago) which has little im- 
pact on the Ha-to-FUV ratios and recent SFHs (< 1 
Gyr ago) which are characterized by a mix of constant 
SF and bursts. We have modeled the older stellar popu- 
lations that dominate a galaxy's total stellar mass using 
a 1 Gyr burst of SF 8 Gyr ago to populate the red gi- 
ant branch and ancient main sequence, as well 1 Gyr of 
constant SFR from 1-2 Gyr ago to account for the con- 
tribution of asymptotic giant branch stars. The precise 
epoch of the ancient and intermediate age SFHs do not 
significantly change any of our analysis. Changing the 
ages of these older bursts results in small variations in 
the absolute values of log(I]_R), i.e., < 0.5 dex. How- 
ever, because the same ancient and intermediate SFHs 
are used for all models, there is no effect on the relative 
values of log(S/j). Finally, to match findings from mea- 
sured SFHs in nearby dw arf galaxies (e.g.. Dolphin et al.l 
120051 iWeisz et al.l 1201 il ). we scaled the SFHs'such that 
SF from > 1 Gyr ago accounts for 90% of the total stellar 
mass formed in each model. 
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Figure 4. An illustrative sampling of select model SFHs we con- 
sider in this paper, plotted over a 500 Myr interval. In panel (a) 
we have denoted the period (P) of a SFH cycle by the grey shaded 
region, the duration (D) of a burst by the red-hatched region, and 
the amplitude (A) is the increase in SFR above the constant equi- 
librium value. Throughout the paper we refer to the models with 
nomenclature such as P125, D40, A30, which corresponds to a pe- 
riod of 125 Myr, a duration of 40 Myr, and amplitude of 30, the 
model SFH in panel (a). Panel (f) is an example of a 'gasping' SFH. 
In to tal we have constructed ~ 1500 model SFHs as described in 

El 
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Figure 5. The model SFHs that best describe the Ha-to-FUV 
flux ratio distribution in order of decreasing bins of stellar mass ((a) 
log(M*) >10.0, (b) 10.0>log(A/*)>9.0, (c) 9.0>log(Af*)>8.0, (d) 
8.0>log(M*)>7.0, (e) log{J\/*) <7.0; see Table Hll, as determined 
by a two sided KS test. The inset shows the evolution of the Ha- 
to-FUV SFR ratio as a function of time for each model SFH. The 
red, green, blue, magenta, and cyan color coding each refer to the 
best fit SFH models in the bins of highest to lowest stellar mass, 
respectively. This color coding scheme will be used throughout the 
remainder of the paper in reference to a particular SFH model. 

We constructed the recent model by assuming a con- 
stant SFR over the past 1 Gyr and superimposing bursts 
of SF. We assume the bursts to be a 1 Myr long episode 
of constant SF, which we then interpolated onto a finer 
time grid to track the evolution of Ha and UV fluxes 
over longer timescales. This basis model was then shifted 
in time, multiplied by the appropriate amplitude, and fi- 
nally linearly co-added with other basis models to obtain 
the desired recent SFH. 

Each model has a specified burst amplitude (A; the 
ratio of maximum SFR to constant SFR), burst dura- 
tion (D), and period between bursts (P). We addition- 
ally characterize our models by the quantity (DxA)/ P, 
which is the ratio of stellar mass formed during a sin- 
gle burst to the mass formed assuming the baseline SFR 
over one period. Selected permutations of A (2, 5, 10, 15, 
20, 30, 40, 50, 75, 100, 150), D (3, 6, 12, 18, 25, 30, 40, 
50, 60, 75, 100, 150, 200 Myr), and P (5, 10, 15, 20, 25, 
30, 40, 50, 60, 75, 100, 125, 150, 175, 200, 225, 250 Myr) 
combined with the constraint P > D resulted in 1441 dif- 
ferent model SFHs. We added 25 more models in which 
we allowed the SFR between bursts to go to zero, i.e., 
the burst amplitude is infinite. We also consider a con- 
stant SFH model (i.e., A = 1). Many of the models we 
consider fall under th e literature parlance of 'bursting' or 
'gasp ing ' SFHs (e.g.. lAnnibah et al.l[200l iMeurer et al.l 
|2009() . Although the definition of both terms is not uni- 
versal, here we consider a 'bursting' model to have ^ > 1 
and P >> D, while for a 'gasping' model D ^ P with 
the deficit in SF equal to 1/A. That is, a gasp is essen- 
tially a decrease in the SFR from an otherwise constant 
rate. We show a select sample of model SFHs in Figure 

H 

Synthesizing the evolution of Ha, FUV, and i?-band 
fluxes was done consistently with the construction of the 
model SFHs. A 1 Myr long burst of SF was input into 



the synthesis code of iBruzual fc Chariot] (|2007f) . The re- 
sulting spectra were convolved with the appropriate filter 
response functions, and the resulting fluxes were linearly 
co-added identically to the SFHs. In calculating these 
integrated quantitie s, we selected the Padova st ellar evo- 
lution models (Bre ssan et al.l[l993l : iFagotto et al. 1994). 
solar metallicity, and a fully populated Chabrier IMF 
with mass limits of 0.1 and 100 M-q. For consistency with 
the observations, we convolved the UV spectrum with 
the GALEX FUV filter response function. The Ha flux 
is modeled as the number of ionizing photos divided by 
10^^ *^, the expectation value for Case B recombination. 
We do not include losses due to the leakage or absorption 
of ionizing photons by dust within the Stromgren sphere. 
In these models, the Ha and FUV fluxes provide for con- 
sistent measures of SF when averaged over the course of 
1 Gyr, i.e., the ratio of the sums, X](Ha)/X](UV), is con- 
stant, even though many individual points may deviate 
from this average. 

Throughout the paper, we assume a fiducial value 
of Ha-to-FUV ratios based on the scenario in which a 
galaxy which has been forming stars at a constant rate 
over its entire lifetime. For clarify in plotting, we have 
set the logarithm of the fiducial Ha-to-FUV ratio equal 
to zero. Consequently, all model and observed Ha-to- 
FUV ratios have been adjusted appropriately such that 
\og[F{H^)lF{FUV)]pioued = \og[F{H^)/F{FUV)] + k, 
where k, = —13.17 

To derive the model i?-band surface brightness values, 
we converted the i?-band fluxes to a surface brightness 
assuming an area of 1 kpc^. In practice, this areal nor- 
malization will increase for more massive galaxies. How- 
ever, as we are only interested in variations in the log- 
arithmic surface brightnesses, altering the areal normal- 
ization by a factor of a few has minimal impact on the 
placement of the models. Overall, this method results 
in model i?-band fluxes (and surface brightnesses) that 
are consistent with observations of galaxies with these 
typical SFRs. 

Uncertainties in the observed Ha, FUV, and _R-band 
fluxes are not inherently factored into the modeled fluxes. 
As this may contribute to the observed scatter, we ap- 
plied independent small random shifts based on observa- 
tional uncertainties in the Ha and FUV fiux es listed in 
iKennicutt et al.l ()2008f ) and lLee et all (|201lD . to each of 
the model fiux values. Specifically, we drew a set of er- 
rors from a randomly generated normal distribution with 
a mean of and a dispersion of 0.05 mag. 

4. RESULTS 

4.1. Comparing Modeled and Observed Ha-to-FUV 
Ratios 

The goal of our analysis is to identify a set of charac- 
teristic SFHs that can match the observed distribution 
of Ha-to-FUV ratios within bins of galaxy stellar masses, 
without violating other observational constraints. To 
determine this characteristic SFH, we compare the ob- 
served distribution of Ha-to-FUV ratios with predicted 
distributions from each of the model SFHs. 

We first divide the observational sample into bins of 
stellar mass, under the assumption that galaxies with 
comparable stellar masses are likely to share a character- 
istic SFH. We chose a simple scheme that consists of five 
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Model KS Test Probabilities 



Ijlll 


iviOQei 


r^erioQ 


Uiira,tion 


AmplituQe 


DxA 
P 


±VD r^rODdDlilLy 


lOg(M0j 




(Myrj 


(Myrj 
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\ ' ) 


log(M*) >10.0 


Best 


125 


3 


10 


0.24 


72.3 




10th Best 


75 


12 


2 


0.32 


65.6 




10th Worst 


125 


12 


150 


14.4 


4.9x10-1° 


10.0>log(M*)>9.0 


Best 


175 


6 


5 


0.17 


99.9 


10th Best 


225 


6 


5 


0.13 


99.6 




10th Worst 


75 


3 


100 


4.00 


2.53x10-12 


9.0>log(M*)>8.0 


Best 


250 


60 


5 


1.20 


94.3 




10th Best 


225 


40 


5 


0.89 


83.4 




10th Worst 


75 


3 


150 


6.00 


1.67x10-31 


8.0>log(M*)>7.0 


Best 


250 


30 


30 


3.60 


97.5 




10th Best 


250 


25 


40 


4.00 


73.4 




10th Worst 


75 


12 


100 


16.0 


1.1x10-22 


log(Af*) <7.0 


Best 


200 


40 


30 


6.00 


99.8 




10th Best 


150 


30 


20 


4.00 


95.8 




10th Worst 


250 


6 


2 


0.05 


3.6x10-^ 



Note. — Parameters and KS probabilities for the best, lOtli best, and lOtli worst fit model SFHs 
in eacli bin of stellar mass. 



bins with sizes of Alog(Af,) ~ 1. The resulting division 
provides a reasonable balance between grouping galax- 
ies of similar masses and populating each bin with an 
adequate number of data points. The number of galax- 
ies, median Ha-to-FUV ratio, and scatter in the ratio for 
each bin are listed in Table [3l 

As a metric for comparison between the models and ob- 
served Ha-to-FUV ratios, we have adopted a two-sided 
Kolmogorov-Smirnov (KS) test. Unlike other statistical 
tests (e.g., i-test, Anderson-Darling), the KS test makes 
no assumption about the form of the parent distribution 
of either population, i.e., they need not be normally dis- 
tributed. Indeed, the distributions of the Ha-to-FUV 
ratios are generally complex and difficult to parameter- 
ize, making the KS test better suited for this particular 
analysis. For simplicity, we will adopt multiples of the 
typical a values to evaluate the quality of the model- 
data agreement (e.g., 1-a ^ 68%, 2-cr ~ 95%, etc). We 
also include parameters and KS test probabilities for the 
10th best and 10th worst fits to highlight the differences 
in the model parameters in the high and low probability 
regimes. 

The results of the model-data comparison are summa- 
rized in Table |4l In each bin of stellar mass, we have 
indicated the highest probability model parameters (i.e., 
P, D, A), the value of {DxA)/P, and the KS probabil- 
ity. In general, in each bin of stellar mass, our highest 
probability models provide reasonable matching to the 
observed distributions of Ha-to-FUV ratios, as judged 
by the KS test. Four of the five bins have models that 
match the data with a probability of > 94%. The lowest 
mass bin {log{M^,) <7.0) and second highest mass bin 
(10.0>log(iV4)>9.0) have KS probabilities of 99.8% and 
99.9%, respectively. The lowest probability model is in 
the highest mass bin (log(M*) >10.0), which has a KS 
probability of 72.3%. Despite the simplicity of the SFH 
models considered, we find all highest probability models 
agree with the data at levels of > 1-3 a. 

Examining the highest probability model parameters 
in each stellar mass bin, we see the SFHs of higher 
mass systems are generally characterized by constant 
SFHs with an occasional modest burst. Specifically, the 



two highest mass bins favor models of short duration 
{D ^ 6 Myr) and relatively modest amplitudes {A < 
10). These SFHs are predominantly constant SFRs, with 
bursts interspersed ~ 5% of the time. The central bin 
{9.0>\og{^I^,)>8.0) is best characterized by a long period 
{P — 250 Myr) of constant SF and an increase in am- 
plitude by a factor of 5 for a 60 Myr duration. The two 
lowest mass bins have highest probability models with 
higher amplitude bursts {A ~ 30), relatively long du- 
rations (D - 30-40 Myr), and large periods {P = 250 
Myr). The highest probability SFH models along with 
the corresponding temporal evolution of the Ha-to-FUV 
ratios are shown in Figure \5\ 

As a rough proxy for the robustness of the highest 
probability model, we rank the models according to the 
probability of matching the data, as judged by the KS 
test. We then examine 10th most probable models in 
each stellar mass bin and find they have similar param- 
eters to the corresponding highest probability models. 
The largest contrast between the best and 10th best 
models is in the highest mass bin, where the duration 
increases by a factor of 4 (from 3 to 12) and the ampli- 
tude decreases by a factor of 5 (from 10 to 2). However, 
the values of {DxA)/P for the best and 10th highest 
probability models are comparable (0.24 and 0.32). This 
simple comparison offers some assurance that we have 
not simply picked out a single model that by chance de- 
scribes the observed distribution. 

Along similar lines, we also consider the 10th low- 
est probability models in each mass bin. These models 
all have vanishingly small probabilities and exhibit pa- 
rameters that are drastically different from the highest 
probability models. For example, the highest probability 
models favor relatively large bursts in the lowest mass 
bin, while the 10th worst model is characterized by long 
periods of constant SF, with bursts of amplitude 2 oc- 
curring only ^ 3% of the time. Conversely, in higher 
mass bins where highest probability models are primar- 
ily characterized by constant SFHs and modest bursts, 
the 10th lowest probability models have typical parame- 
ters of short duration, high amplitude bursts. 

An examination of full parameter space reveals a de- 
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Table 5 

Model and Observed Ha-to-FUV Ratios Relative to the 
Fidueial 



Bin 


Above 


Above 


Below 


Below 




Observed 


Model 


Observed 


Model 


(1) 


(2) 


(3) 


(4) 


(5) 


log{A/*) >10.0 


31.3 


30.6 


68.7 


69.4 


10.0>log{A/*)>9.0 


35.0 


36.5 


65.0 


63.5 


9.0>log(M*)>8.0 


29.1 


31.5 


70.9 


68.5 


8.0>log(M*)>7.0 


15.5 


13.6 


84.5 


86.4 


log{J\/*) <7.0 


21.7 


16.3 


78.3 


83.7 



Note. — A comparison between the observed and pre- 
dicted fraction of galaxies above and below the fiducial, 
per stellar mass bin. Columns (2) and (4) correspond to 
the observations while Columns (3) and (5) are values for 
the models. 

gree of degeneracy in model parameters. Specifically, 
models with similar values of {DxA)/P also have com- 
parable KS probabilities. Thus, we consider the highest 
probability model as a representation of a set of models 
with similar values of (DxA)/ P. 

4.2. Evolution of Predicted Fluxes, Masses, and SFRs 
4.2.1. Model SFRs versus Stellar Mass and Sjj 

As demonstrated inl Meurer et al.l (|2009D . a model SFH 
that matches the distribution of Ha-to-FUV ratios may 
not be able to simultaneously account for other observ- 
ables, such as i?-band surface brightness or stellar mass. 
With this in mind, in Figure IHl we compare the evolution 
of the Ha-to-FUV ratios versus stellar mass (top) and 
i?-band surface brightness (bottom) over one burst cy- 
cle against the observed data. For illustrative purposes, 
we have excluded the simulated observational noise from 
this particular figure. 

The relative evolution of stellar mass and surface 
brightness have been derived entirely from the models. 
We centered the starting point of the models on the me- 
dian values of the observations. This choice resulted 
in overlap between points from different models, and in 
those cases, small adjustments (< 0.3 dex) in the ab- 
solute placement of mass or luminosity were made to 
improve clarity. Loops are composed of discrete model 
points sampled at time intervals of 0.1 Myr. For ref- 
erence, we have added time labels that correspond to 
distinct events in the burst cycle to the loop in the low- 
est surface brightness regime in Figure [6] (cyan points), 
which has a period of 250 Myr, a burst duration of 40 
Myr, and an amplitude of 30 (P250, D40, A30). 

During a burst, the Ha and FUV SFRs are not in equi- 
librium. Early in the burst, the production of Ha pho- 
tons exceeds that of the FUV. However, stars that are 
FUV bright, are less massive, and have longer typical 
lifetimes than those responsible for Ha production. The 
difference in stellar lifetimes leads to a build up of FUV 
emission over the course of the burst, resulting in a net 
decline in the Ha-to-FUV ratio. We see this behavior 
concretely in the magenta 'loop' in Figure |6l Immedi- 
ately following the beginning of a burst, (denoted by '0 
Myr'), the Ha-to-FUV ratio is at its peak, as expected. 
As the burst continues over the subsequent duration of 
40 Myr, the Ha-to-FUV ratio steadily declines. 

Following the end of the burst, we see a significant de- 
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Figure 6. The predicted Ha-to-FUV ratios from one burst cycle 
of the best fit SFHs plotted over the observational data. The red, 
green, blue, magenta, and cyan color coding each refer to the best fit 
SFH models in the bins of highest to lowest stellar mass, as detailed 
in Figure [5] To help illustrate the differences in the fiux ratios, a 
constant k has been added. The value of n is —13.17, the negative 
of the expected fiux ratios from our models for a constant SFH. To 
illustrate the evolution of the different quantities, we have plotted the 
models without simulated observational noise. The relative evolution 
in stellar mass and surface brightness of the loops has been determined 
entirely by the model SFHs. In the bottom panel, the age labeling 
(black) corresponds to time after the burst cycle has started. The 
model points arc plotted for equal time intervals of 0.1 Myr. The 
horizontal black dot-dashed line in each panel is the fiducial ratio. 



crease in the Ha-to-FUV ratio between 40 and 50 Myr. 
There is negligible Ha flux after the burst, as the more 
massive stars expire within a few Myr of the termination 
of star formation. Thus, the build up of FUV emission 
drives the net Ha-to-FUV ratio to lower values. In ad- 
dition, the SFH has returned to a constant, lower level 
SFR during this 10 Myr period, leading to an equilibrium 
state for the production of Ha and FUV fiux. 

Finally, from 50 to 250 Myr, we see a steady increase in 
the Ha-to-FUV ratio. The FUV excess, built up during 
the burst, continues to fade, and the contribution of the 
constant SFR to the net Ha-to-FUV ratio increases. At 
250 Myr, near the end of the burst cycle, the Ha-to-FUV 
ratio approaches the equilibrium value. 

We now examine the same sequence, only for the evo- 
lution of i?-band surface brightness. During the burst, 
the i?-band surface brightness increases due to the pro- 
duction of luminous young stars, reaching a maximum at 
the end of the burst. Relative to the initial value, the sur- 
face brightness has increased by a value of log(E/j) ~ 0.5. 
For the subsequent 10 Myr, the _R-band surface bright- 
ness decreases. The death rate of luminous stars exceeds 
the birth rate, meaning the galaxy cannot sustain the 
maximum surface brightness level obtained during the 
burst phase, and consequently begins to fade. Finally, 
from 50 to 250 Myr, the i?-band surface brightness con- 
tinues to fade and approaches the value consistent with 
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Figure 7. The same as Figure [6l only including the efTects of 
simulated observational noise. 
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Figure 8. The same as Figure [71 only now plotted versus 
SFR(Ha) and Mb- 



that expected from a constant SFH. 

4.2.2. Flux Evolution with Simulated Noise 

In Figure [71 we now consider the model loops including 
the effects of simulated noise, introduced to mimic un- 
certainties in the observed fluxes. The inclusion of noise 
adds to the dispersion in the models points. The model 
points represent equal time intervals of 0.1 Myr, such 
that the density of points is directly proportional to the 
amount of time spent at a given phase of the burst cycle. 

From this information, we conduct a simple compari- 
son between the fraction of model and data points above 
and below the fiducial Ha-to-FUV ratio. As shown in Ta- 
ble [H the predictions from the models agree with the ob- 
servations within 5%. In general, ~ 30-35% of the ob- 
served massive galaxy sample is above the fiducial. The 
corresponding models predict a similar fraction. Obser- 
vations of lower mass systems show that ^ 15-20% of the 
galaxies are above the fiducial. The highest probability 
SFH models in these bins also occupy a similar range, 
indicating that models of galaxies bursty SFHs naturally 
account for low mass galaxies with both low and normal 
Ha-to-FUV ratios. 

4.3. Ha-to-FUV Evolution vs. SFR(Ha) and Mb 

In Figure [HI we plot the Ha-to-FUV ratios versus 
SFR(Hq!) (top) and versus Mb (bottom), including the 
effects of simulated noise as in >i4.2.2l In the top panel, 
we see that the predicted change in the Ha-to-FUV ratio 
shows an expected strong correlation with the Ha based 
SFR, such that galaxies with higher Ha-to-FUV ratios 
generally have higher Ha SFRs. From the observations, 
we see that galaxies with low Ha SFRs and low Ha-to- 
FUV ratios are quite rare, which agrees with the SFH 
model predictions. 

Plotted versus Mb, the model SFHs indicate an in- 
crease of ~ 2 mag in blue luminosity during a burst. The 



majority of the time, however, the model SFHs show that 
galaxies are in lower luminosity states with Ha-to-FUV 
ratios below the fiducial. 

5. DISCUSSION 

5.1. Comparison with Previous Studies of Ha-to-FUV 

Ratios 

Our study is not the first to model the effects of SFHs 
on Ha-to-FUV ratios. Within the past decade three no- 
table studies (Ug lesias- Paramo et all 12004 iBoselli et al.l 
120091 iMeurer et' al.li2009t ) have compared Ha and FUV 
fluxes from sets of SFH models with those from ob- 
servations. These three studies represented a forward 
step in considering a range of SFH parameter space 
beyond past analysis, which were generally restricted 
to eith er constant or instantaneous burs t models (e.g. , 
Buatet al. 1987; Buat 1992; Glazcbro ok et al.l 119991: 
Yan ct al. 1999; Sullivan ct aL| l2000t iBell fc KennicuttI 
2001; Moorwood et al. 2000). 

Likewise, the analysis presented in this paper pro- 
vides another incremental step toward understanding 
the impact of non-constant SFHs on observed Ha-to- 
FUV ratios. The primary differences between this study 
and th ose o f [iglesias-Par a mo et al.. (2004), Boselli et al] 
()2009f ). and IMeurer et al.! (f2009l ) are in the construction 
and analysis of model SFHs and in the method of cor- 
recting for internal dust attenuation. In this section, we 
will focus on the differences in model SFHs. 

The model SFHs considered in this paper represent a 
fine sampling of SFH parameter space. In contrast, two 
of the three previous studies consider a coarse sampling, 
genera lly intended to brac ket ranges of SFHs. For ex- 
ample, IMeurer et al.l ()2009D consider 18 single burst and 
gasp model SFHs with durations of 10, 10, 1000 Myr 
and amplitudes of 2, 10, and 100, and find that none are 
able to simultaneously account for the observed trends of 
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Table 6 

Model and Observed Ha-to-FUV Ratios Relative to the 
Fidueial 



Bin 


Above 


Above 


Below 


Below 




Observed 


Model 


Observed 


Model 


(1) 


(2) 


(3) 


(4) 


(5) 


log{A/*) >10.0 


31.3 


30.6 


68.7 


69.4 


10.0>log{A/*)>9.0 


35.0 


36.5 


65.0 


63.5 


9.0>log(M*)>8.0 


29.1 


31.5 


70.9 


68.5 


8.0>log(M*)>7.0 


15.5 


13.6 


84.5 


86.4 


log{J\/*) <7.0 


21.7 


16.3 


78.3 


83.7 



Note. — A comparison between the observed and pre- 
dicted fraction of galaxies above and below the fiducial, 
per stellar mass bin. Columns (2) and (4) correspond to 
the observations while Columns (3) and (5) are values for 
the models. 



Hg-t o-FUV ratios versus Ha and_/?-band surface bright- 
ness. Ilglesias-Paramo et al.l ()2004[ ) considered a slightly 
broader range SFH parameters and analyze 36 model 
SFHs, but reach an ambiguous conclusion concerning the 
influence of SFHs on Ha-to-FUV ratios. Our SFH mod- 
eling finds that none of the extreme models, i.e., those 
considered bracketing the range of parameter space, pro- 
vide good descriptions of the data. Thus, it is not en- 
tirely surprising that these prior studies did not find well 
matched model SFHs. 

The 'micro-SFHs' considered by IBoselh et al.1 (pOOl ) 
represent a limited, but fine sampling of SFH parameter 
space. The focus of micro-SFHs is on the effects of the 
more recent burst on the Ha luminosity. Episodic bursts 
with durations of 2, 5, and 10 Myr over the course of 
100 Myr were considered to simulate the effects of vari- 
ous HIT regions turning on and off. The principle behind 
this method is similar to our analysis, except we consider 
a larger number of possible model permutations and the 
contribution of ancient stellar populations, allowing us 
to accurately assess the impact of recent SFHs on stel- 
lar mass and surface brightness, as well as the observed 
Ha and FUV fluxes. We further note that the models 
considered bv lBoselli et al.l (|2009t ) are similar to those in 
our highest mass bins, but they do not consider mod- 
els with parameters such as those in our lowest mass 
bins. This se ems appropriate, however, as the sample of 
IBoselli et alj ([2009) does not extend to comparably low 
luminosity galaxies. 

A further point of contrast between our study and the 
previous three is the division of the sample for analysis. 
Each of the previous studies attempt to find a single SFH 
that best describes the entire observed sample. By bin- 
ning the data in increments of stellar mass, we do not 
impose the restriction that a single model SFH must be 
able to explain the observed Ha-to-FUV ratios in both 
massive galaxies such as M51 and low mass dwarfs such 
as GR8. Indeed, as we have shown, the highest probabil- 
ity models converge on drastically different characteristic 
SFHs for galaxies of different masses. 

5.2. SFHs and the IMF 

An underlying assumption to any calibration of SFR 
indicators is the nature of the IMF. In the conversion 
from observed flux to SFR, the IMF is assumed to be 
universal, i.e., the same with respect to time and envi- 
ronment, and fully populated, i.e., the relative number 



of stars per unit mass interval is always the same. How- 
ever, whether either of these criteria are true in nature 
remains a he avily debated an d open question (see the 
review by Bas tian et 311120101 ). 

Unfortunately, only a handful of variable IMF scenar- 
ios in the literature have been discussed in the context 
of Ha-to-FUV ratios. Extracting results from some of 
these studies, we proceed to discuss these findings along- 
side predictions from the model SFHs analyzed in this 
paper. 

5.2.1. Stochastically Sampling the IMF 

We first consider the case of a universal, but not fully 
sampled IMF. The effect of this stochastically sampled 
IMF sce nario on the Ha- to-FUV ratio has been simu- 
lated in iLee et al.l ()2009b[ ). Assuming a constant SFH, 
the authors found that the predicted turnover in Ha-to- 
FUV ratios occurred at too l ow an Ha luminos^ ity relative 
to the data (see Figure 7 in ILee et al.ll2009bl) . That is, 
randomly sampling a universal IMF does not appear to 
adequately account for the observations. 

This res ult is c onfirm ed by more recent simulations 
in da Sil va et all (|2011| ). |Fumagalh et al. (2011), and 
lEldridga (|2011h . In each paper, the authors verify the in- 
ability of random sampling of the IMF to account for the 
observed trend in Ha-to-FUV ratios, under the assump- 
tion of a constant SFH and a fully populated cluster mass 
function. Instead, each study demonstrates that stochas- 
tic sampling of both the stellar and cluster mass func- 
tions are necessary to produce Ha-to-FUV ratios that 
are consistent with observations. The effect of stochas- 
tically sampling the cluster mass function is that lower 
mass galaxies form massive clusters less frequently than 
higher mass systems. Along the same lines, stochasti- 
cally sampling stellar IMF implies that high mass stars 
have a higher probability of forming in high mass clus- 
ters. In light of these scenarios, we would then expect 
low mass galaxies to form fewer high mass clusters, and 
hence fewer high mass stars when compared to more mas- 
sive galaxies. Therefore, low mass systems would have 
fewer massive stars capable of producing Ha, and the 
Ha-to-FUV ratios in low mass systems is expected to 
be less than for more massive galaxies. Crucially, these 
models posit that the formation of massive clusters in 
low mass galaxies are rare, but not impossible, which is 
the case for certain prescriptions of a variable IMF (see 
g5221)- More detailed discussion and quan titative analy - 
sis of these effec t s can be fo und in da Sil va et al.l (|2011| ). 
iFumagalli et~all (|20Tll) . andHdridga (201llV 

5.2.2. SFHs and the Integrated Galactic IMF 

The second possibility we consider is that the IMF is 
not universal. This scenario has been widely debated 
in the literature, but a co nsensus has yet to be reached 
(e.g.. iBastian et al.ll20Tot) . We focus on one particular 
scenario , the Integrated Ga lactic IMF (IGIMF) as de- 
tailed in lWeidner fc Kroupal (2005), which makes specific 
predictio ns for Ha-to-FUV ratios as functions of Ha lu- 
minosity (|Pflamm-Altenburg et~a l. 2007, 2009). 

There are three tenets of the IGIMF: (i) all stars are 
formed in clusters; (ii) the maximum mass of a cluster 
star is determined by the mass of the cluster in which it 
is formed; (iii) the maximum cluster mass is a determin- 
istic function of the integrated SFR of a galaxy. From 
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these three postulates, the IGIMF predicts a systematic 
increase in the ratio of Ha-to-FUV fluxes as a function 
of Ha luminosity, under the as sumption of a cons tant 
underlying SFH (Pflamm-Alten burg et "alll2007l [2009). 

In what follows, we will compare predictions from the 
IGIMF and model SFHs with our observations. First, 
we will consider the standard case of the Ha-to-FUV 
ratios as a function of Ha luminosity. The IGIMF pre- 
dictions in this case are t aken d irectly from the work of 
IPflamm -Alten burg et all (|2009D . We will then compare 
the model predictions to the data as a function of galaxy 
stellar mass. In this case, to compute the galaxy stellar 
mass for the IGIMF models, we have assumed a con- 
stant SFH over the history of the galaxy at a level of 
the current Ha SFR. We discuss on the validity of this 
assumption below. 

As a firs t compari s on, we follow the discussion of the 
IGIMF in iLee et all (j2009b[ ). and over-plot predictions 
for the 'standard' and 'minimuml' IGIMF models (solid 
and dashed navy lines) on the observed Ha-to-FUV ra- 
tios versus Ha luminosity in the top panel of Figure [9l 
The two IGIMF models correspond to different convolu- 
tions of the cluster mass function with a stellar IMF, such 
that the minimall model is the least deviant from a stan- 
dard Salpeter IMF (see IPflamm-Altenburg et all 120071 
120091 for speci fic details). We have adopted t he IGIMF 
models from IPflamm-Altenburg et al.l (|2009D that as- 
sume a constant SFH over the history of the universe. 

We have also over-plotted the predictions from the 
highest probability model SFHs. Each of the colored 
points represents the median Ha-to-FUV ratio from 
models shown in Figure [5] and listed in Table S] The 
placement on the horizontal axis corresponds to the me- 
dian Ha luminosity in each mass bin for the observed 
sample. The error bars represent the 5th and 95th per- 
centiles (~ 2-cr) in the model distributions. 

In the top panel of Figure IHl we see that both the 
IGIMF and SFH models qualitatively agree with obser- 
vations. At high Ha luminosities, both types of model 
indicate that ratios should be near the fiducial. For de- 
creasing Ha luminosity, both models predict a decline 
in the Ha-to-FUV ratio that is in reasonable agreement 
with the mean observational trends. Of the two IGIMF 
scenarios, the 'minimall' model appears to more closely 
follow the mean of the data better than the 'standardl' 
model, particularly when the data are considered as a 
function of stellar mass. 

The scatter in the Ha-to-FUV ratios is another quan- 
tity we consider. For Ha luminosities > lO'^*'^, we see 
a number of data points that lie within the range of the 
SFH model predictions, but above the two IGIMF curves. 
We note that it is not clear how the IGIMF can explain 
galaxies that lie significantly above the predicted curves, 
and pursue a more detailed discussion of this point below. 

In the bottom panel of Figure[3J we repeat the compar- 
ison outlined above, only this time we plot galaxy stellar 
mass on the x-axis. Again, we see that both the SFH 
and IGIMF models predict a general decline in the Ha- 
to-FUV ratios as a function of decreasing stellar mass. 
However, both IGIMF models tend to under predict the 
Ha-to-FUV ratios for galaxies < 10* M©. In comparison, 
the median values of the SFH models do not decrease as 
sharply as those of the IGIMF, and instead exhibit a 
larger dynamic range of possible model values. 
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Figure 9. The observed Ha-to-FUV flux ratios {grey points) versus 
Hq luminosity (top panel) and stellar mass (bottom panel). The large 
eolorcd point correspond to the median values from the best fit SFH 
modes as described in i|4.1l The error bars on the models correspond 
to the 5th and 95th percentiles of the Ha-to-FUV flux distributions for 
each model SFH. The dashed and solid navy lines correspond to pre- 
dictions of the IGIMF 'minimall' and 's t andard' mo dels, respectively, 
as described in IPflamm-Altenburg et al.l 1 120071 . 120091) . The black dot- 
dashed line represents the fiducial ratio for our models. Axes repre- 
sented Ha and FUV SFRs are grey as they are not valid in regions of 
low SFR(Ha) according to the IGIMF models. The two model types 
are indistinguishable in the upper plot when the flux ratios is consid- 
ered versus Ha luminosity. However, when plotted versus galaxy stellar 
mass, the IGIMF models are no longer consistent with the majority of 
points at low masses, which lie above the model curves. 



Of particular interest is the relationship between the 
two types of models and the scatter in the data. The 
broad dispersion of Ha-to-FUV ratios predicted by the 
model SFHs are generally in good agreement with the 
observations. The most notable discrepancies are the 
handful of galaxies that have Ha-to-FUV ratios that are 
lower than predicted by the models. This may be an 
indication that stochastic IMF effects need to be included 
to explain the extremely low Ha-to-FUV points. 

The IGIMF models are generally not in good agree- 
ment with the observed scatter. In the best case, the 
IGIMF models track the decline of average Ha-to-FUV 
ratios as a function of decreasing Ha luminosity, but 
there are a number of galaxies with ratios significantly 
below the IGIMF predicted curves, and a smaller frac- 
tion that lie above. The disagreement worsens when the 
comparison is done as a function of stellar mass. In par- 
ticular, for galaxies with stellar masses < 10* M©, there 
are a substantial number of galaxies that have Ha-to- 
FUV ratios significantly above the IGIMF predictions. 

These data points are challenging to understand in the 
framework of the IGIMF. Take n at face value, the IGIMF 
is a d eterministic theory (e.g., IPflamm-Altenburg et al.l 
l2009f ). That is, it makes specific predictions about the 
the functional form of the Ha-to-FUV ratio, namely that 
Ha should be deficient in galaxies with low masses or low 
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SFRs. Yet, the data show a number of low mass systems 
with normal or high Ha-to-FUV ratios, implying the fre- 
quent presence of massive star(s) in low mass systems; 
a observational finding that does not seem compatible 
with the tenets of the IGIMF. 

Admittedly, there are two areas of uncertainty with 
this interpretation. First, the connection between to- 
tal galaxy stellar mass and the IGIMF is not well de- 
fined in the literature. Although such a connection 
has been invoked in other analyses of the IGIMF (e.g., 
IWeidner & Kroupa 2005; Elmegreen 2006; Bosclli"eFaI] 
[2009), extrapolating the current Ha SFR to a constant 
lifetime SFH to determine the stellar mass is not a valid 
assumption. In the event that we have over-estimated 
the galaxy stellar mass, the IGIMF curves would shift 
to the left in the bottom panel of Figure O In order to 
match the mean observed flux ratios, the extrapolated 
galaxy masses need to be decreased by a factor of ^ 10. 
Even with this correction applied, the IGIMF curves still 
do not match the observed scatter in the data very well. 

The second area of uncertainty is the connection be- 
tween the IGIMF and assumed form of the SFH. Both the 
'standard' and ' minimal 1' IGI MF models assume acon- 
stant SFH (e.g., iPflamm-Altenburg et al.l [20071 [2001 ). 
Including the effects of a bursty SFH may be able to pro- 
duce Ha-to-FUV ratios higher than currently predicted. 
As seen from the model SFHs, during the early phases 
of a burst, the Ha-to-FUV ratios are near their peak. 
Combining this mechanism with the IGIMF provides at 
least one possible explanation for galaxies with higher 
than predicted Ha-to-FUV ratios. 

However, the convolution of bursts with the IGIMF is 
not without complications. First, if bursts are required 
to explain Ha-to-FUV ratios above the IGIMF predic- 
tions, this implies that all such galaxies are currently in 
the midst of a burst. However, it seems unlikely that 
the majority of low mass galaxies in our sample should 
be undergoing simultaneous bursts. Second, under the 
auspices of the IGIMF, one also needs to consider the 
absolute SFR, not just the relative burst amplitudes. At 
low absolute SFRs such as lO""* Mq yr^^, the IGIMF 
model predict that log[F{Ha)/ F{FUV)] is 1. Ex- 
trapolating this SFR over 13.7 Gyr yields a stellar mass 
of ~ 10^ Mq. At this stellar mass, we see a handful of 
galaxies with \og[F{Ha)/F{FUV)] > 0. The only way 
for such low mass galaxies to have these high Ha-to-FUV 
ratios is through intense bursts of SF that significantly 
increase the amount of Ha emission. Additionally, we 
note that such bursts should not be too long in duration 
or too high in frequency or the final mass of the galaxy 
could be substantially increased. 

Overall, we find that bursty SFH models can qualita- 
tively account for the general decline and scatter in the 
observations as a function of decreasing Ha luminosity 
or galaxy stellar mass, whereas the simplest form of the 
IGIMF models do not. We emphasize that this finding 
does not invalidate the IGIMF models, nor suggest that 
a variable IMF is not possible, only that there appear 
to be other viable explanations for the observed Ha-to- 
FUV ratios that do not involve a systematically varying 
IMF. Indeed, we believe that future model efforts that 
include more complex SFHs, and reasonable variations 
on the IMF are essential for understanding SF processes 
in galaxies. We suggest some possible paths forward in 



the following section. 

5.3. Toward More Realistic SFH Models 

In this paper, we have found a set of SFH models that 
match the trend and scatter in observed Ha-to FUV ra- 
tios. However, we caution that the broad applicability of 
our conclusions should be tempered by the inherent sim- 
plicity of the models. SFHs in nature are undoubtably 
more complex than periodic bursts with fixed amplitudes 
and durations. We have considered characteristic SFH 
models for galaxies in comparable mass ranges, but it is 
likely that even galaxies with similar stellar masses do 
not share identical SFHs. Indeed, studies within the Lo- 
cal Group and nearby universe have shown that galaxies 
with similar m asses exhi bit diversity i n their SFHs (e.g., 
IDolphin et an r2005: Wei sz et al.|[20lT[ ). 

Determining the degree to which any model SFH truly 
refiects those in nature is a challenging and presently un- 
resolved question. A comprehensive assessment of the re- 
alistic nature of the model SFHs used in this paper is be- 
yond the s cope of this study, but is the subject of future 
work (e. g., IB. Johnson et aLllin prep.l : ID. R. Weisz et al] 
lin prep.r ). Instead, in this section we outline the limi- 
tations of the simple model SFHs and discuss ways in 
which future SFH model efforts can be improved. 

One increasingly common method of directly measur- 
ing SFHs is from analysis of color-magnitude diagrams 
(C MDs) of resolved stell ar populations in nearby galaxies 
(see lTolstoY et all 12009 and references therein). In Fig- 
ure IIOI we show an example of six such SFHs measured 
from CMDs based on Hubble Space Tele scope imag- 
ing o f dwarf galaxies in the M81 Group (iWeisz et al.l 
mm iDalcanton eFall [20091 : IWeisz et all 1201 ID . The 
solid black and dot-dashed magenta lines represent the 
same SFH binned to different time resolutions, which we 
discuss below. 

These measured SFHs provide an indication of the 
complexity of galaxy-wide star formation. The most 
striking characteristic in these measured SFHs is the de- 
gree of stochasticity in the amplitude and duration of 
star formation episodes. Further, there is no clearly de- 
fined period associated with any of the measured SFHs. 
There is also a conspicuous lack of stochasticity in the 
model SFHs (e.g.. Figure [J). Clearly, non- uniform pa- 
rameterization of bursty SFHs is an important next step 
toward constructing more realistic SFHs. 

Of course, there is some degree of similarity between 
the modeled and measured SFHs. Focusing on the most 
recent 100 Myr (insets in Figure [TUjl . where the time res- 
olution of the CMD-based SFHs is the highest, a hand- 
ful of bursts with amplitudes of ~ 10-50 are present in 
the measured SFHs, similar to those found in the high- 
est probability models for low mass galaxies. The mea- 
sured SFHs also suggest that extended duration, high 
amplitude bursts (e.g., D > 100 Myr, A > 50) are 
probably not common. Burst durations of hundreds of 
Myrs have bee n measured in a mino rity of starburst 
dwarf galaxies (jMcQuinn et al.l I2010ai rb[) , but with av- 
erage burst amplitudes of less than a factor of 5. How- 
ever, these global, long duration bursts are superimposed 
with a short timescale 'flickering' of enhanced SF, which 
result in short-term burst amplitudes > 10 on 10 Myr 
timescales. 

Average burst amplitudes in these systems are typi- 
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Figure 1 0. Select SFHs m oasurcdf from analysis of HST-based opti- 
cal CMDs IWeisz et al.l200a) . The SFHs plotted as blaek solid lines are 
shown at the highest possible time resolution, while the dot-dashed ma- 
genta lines indicate the more conservative time binning scheme adopted 
in Wcisz ct al. (2008[). The insets show the most recent 100 Myr of each 
SFH, where the CMDs provide for high time resolution on the derived 
SFHs. Errors on the SFRs have been omitted for clarity, but the in- 
crease in time bin size generally corresponds to smaller uncertainties 
in the SFRs. We draw attention to two features: the amplitudes of SF 
events and the complexity of the measured SFHs. Characteristic am- 
plitudes in the measured SFHs are qualitatively similar to those in the 
best fit model SFHs, but clearly a more precise analysis is needed to 
quantify the similarity. Relative to the simple model SFHs we consider, 
the measured SFHs show stochastic variations in both the frequency 
and duration of SF events, as well as their amplitudes. Incorporat- 
ing such features into future SFH models in an important step toward 
making them more realistic. 



cally a factor of a < 5 over the course of hundreds of 
Myr, ahhough variations in the SFRs on timescales of 
order ~ 10 Myr can vary by factors of 10 or more. De- 
terminations of exact amplitude and duration parame- 
ters directly from CMD-based SFHs are challenging due 
to the logarithmic time resolution of the measured SFHs. 
We explicitly illustrate the issue at hand in Appendix 1X1 
Integrated fluxes provide further empirical constraints 
for SFH modeling efforts. A realistic set of model SFHs 
should be able to reproduce such observations as the dis- 
tribution of observed Ha EWs, the UV luminosity func- 
tion, and optical broadband colors and luminosities. The 
simple model SFHs we consider are not in strong agree- 
ment with either observed distrib utions of ei t her Ha EW 
or UV luminosit y, as presented in iLee et ahl ()2009af ) and 
iLee et all ()201lD . While both the model and observed 
distributions are well described by lognormal functions, 
the functional form of the model predictions has too 
many galaxies near the median value and not enough 
at the extremes. Some of the difference can be recon- 
ciled by increasing the simulated noise in the models, but 
we would prefer to avoid this type of fine-tuning. Fur- 
ther, the addition of stochasticity would also mitigate 
the over-density of galaxies near the median of the dis- 
tribution. Quantitatively testing these effects is beyond 



the scope of this paper, but is discusse d in the context 
of the SLUG population s ynthesis code ()Fumagalli et al.l 
IMltldaSTlva et al.l[20ll . 

6. SUMMARY 

We have compared predictions from simple, periodic 
SFHs with observed Ha-to-FUV SFR ratios observed as 
part of the IIHUGS and Spitzer LVL programs. From 
a suite of ^ 1500 simple model SFHs, we have identi- 
fied a set of models that matches the observed distribu- 
tion of Ha-to-FUV SFR ratios over a range of galaxy 
stellar masses. We find that high mass galaxies tend 
to have characteristic SFHs that are predominantly con- 
stant, with relatively modest amplitude bursts {A < 10) 
spaced by ~ 100 Myr. We find that lower mass galaxies 
are best described by SFHs with burst amplitudes of ~ 
30, and inter-burst spacings of ^ 150-200 Myr. These 
model SFHs have burst amplitudes that are compara- 
ble to those inferred from CMDs of nearby star forming 
dwarf galaxies. 

Using results from the highest probability models, we 
compared the predicted evolution of the Ha-to-FUV 
ratios versus stellar mass, i?-band surface brightness, 
SFR(Ha), and Mb, and found that in all cases, the mod- 
els were in good agreement with observations. 

We find that highest probability models are well 
matched to the observed systematic decline in the Ha- 
to-FUV ratios as a function of both decreasing Ha lumi- 
nosity and decreasing galaxy stellar mass. 

Variations in the high mass stellar IMF have been pro- 
posed as an explanation for this trend. We therefore 
compare predictions from the model SFHs with those 
from the IGIMF theory, which also makes specific predic- 
tions for Ha-to-FUV ratios. From this comparison we see 
that stochastic sampling of the IMF cannot re-produce 
the observed trend. A convolution of stochastically sam- 
pled stellar and cluster mass functions are necessary to 
match observations. We find that both the SFH and 
IGIMF predictions are in good agreement with the Ha- 
to-FUV ratios as a function of Ha luminosity. However, 
when considered versus stellar mass, it is not clear if the 
IGIMF allow for Ha-to-FUV ratios higher than the pre- 
dicted curves, whereas the model SFHs naturally account 
for the observed scatter. 

Instead, the we find that only time variable SFHs are 
able to explain both the mean trend in the data, as well 
as the scatter. While we have only considered toy model 
SFHs, they share characteristics with the stochastically 
sampled cluster and ste l lar m as s function scena r io pre - 
sent ed inlda Silva et al.l ()2011[ ). iFumagalh eFall (|201lD . 
and|Eldridgd(|2011(). Essentially, at low SFRs, a stochas- 
tically sampled cluster mass function rarely produces 
massive clusters. However, when a massive cluster does 
form, the SFR of the galaxy increases, and can effec- 
tively look like a burst of SFH. While the frequency of 
massive cluster formation in low mass galaxy is not well 
understood, the new population synthesis code, SLUG 
(|da Silva et al.l 120111: |p\imagalh et all 1201 It ) wih permit 
investigations on the impact of stochastic cluster forma- 
tion on the galaxy wide SFH and broadband flux prop- 
erties. 

Finally, we discuss the limitations of our toy model 
SFHs, including the assumptions of periodicity and lack 
of stochasticity. In addition, we find that one major 
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challenge is understanding differences between our model 
predictions and observations of the Ha and UV lumi- 
nosity functions, as well as precisely matching observed 
optical and NIR galaxy luminosities. Exploring these dif- 
ferences are one of the topics that wil l be explo red in the 
upcoming paper of IB. Johnson et"aLl (|in prep.f l. 
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APPENDIX 

A. THE EFFECTS OF LOGARITHMIC TIME RESOLUTION OF CMD-BASED SFHS FOR DERIVING MODEL SFH 

PARAMETERS 

In principle, one can synthesize Ha-to-FUV ratios from CMD based SFHs, eliminating the need for the intermediate 
step of modeling SFHs. However, the non-uniformly linear time resolution of CMD based SFHs presents a challenge to 
such efforts. Much like stellar isochrones, CMD based SFHs are naturally sampled in uniform logarithmic time steps, 
which are not equally spaced in linear time. The combined coarseness of the time binning for lookback times 100 
Myr and the sensitivity of the Ha-to-FUV ratios to the input parameters limit the utility of CMD based SFHs for this 
particular purpose. 

As a simple demonstration of the effects of time resolution , in FigurelTUl we sh ow the same SFHs binned at the highest 
time resolution (solid black) and the resolution adopted by iWeisz et al.l (|200^ (dot-dashed magenta) . Although these 
represent the same underlying SFH, the choice in time-binning gives a strikingly different impression of the galaxy 
SFHs. Extending this exercise to the models considered in this paper, in Figure [TT] we plot the example model SFHs 
from Figure m convolved to the two CMD based time binning schemes. In most cases, the input SFHs are sufficiently 
degraded to an unrecognizable state. There is, of course, some dependence on the particular SFH, as well as the quality 
of the observed CMD. A more detailed discussion of time resolution and SFH recovery from CMD can be found in 
iMcQuinn et"all (|2Q10aD . 
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Figure 11. The sample model SFHs presented in FigurelH resample d at the high e st tim e resolution binning available for CMD 
based SFHs (solid black lines) and the broader binning scheme of IWeisz et al.l l|2008t ) in the magenta dot-dashed line. For 
reference, the original model SFHs have been plotted in the navy dashed lines. The native logarithmic time binning scheme 
for CMD based SFHs can make it difficult to differentiate between various SFH models for times > 100 Myr. Bursts that are 
closely spaced, low in amplitude, or short in duration tend to get washed out and appear as relatively small deviations above 
an otherwise seemingly constant SFH. In addition to challenges posed by the logarithmic nature of CMD-based SFHs, there are 
degeneracies in certain stellar features, e.g., the luminous MS, a nd observational effects such as differential extinction which can 
further reduce the time resolution. See IMcQuinn et al.l lj2010a[ ) for further discussion of CMD time resolution effects. 



